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Abstract 

The definition of spacings associated to a sequence of random variables is 
extended to the case of random vectors in [0, l] 2 . Beirlant & al. (1991) give 
an alternative proof of the Le Cam (1958) theorem concerning asymptotic 
normality of additive functions of uniform spacings in [0,1]. I adapt their 
technique to the two-dimensional case, leading the way to new directions in 
the domain of Complete Spatial Randomness (CSR) testing. 



1 Introduction 

Testing the uniformity of real random variables (r.v.) can be done in sev- 
eral ways : using Chi-square tests, tests based on the empirical distribution 
function (e.d.f.), tests based on spacings ... The latter ones have been exten- 
sively studied (Pyke, 1965) and recommended, for example, for the analysis 
of the local renewal structure of a point process (Deheuvels, 1983a). 

In higher dimensions, when dealing with a spatial point pattern U G S C 
R d , one first wishes to know whether it satisfies the CSR hypothesis : is the 
spatial process governing U a homogeneous Poisson process ? This question 
is equivalent to the following : given the number of points in the pattern (also 
called events), are these points uniformly and independently distributed in 
S (Moller & Waagepetersen, 2004) ? 

I concentrate here on point patterns distributed in rectangles in R 2 , which 
is similar, after linear transformation of the coordinates, to testing the uni- 
formity in [0, l] 2 . 
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Most of two-dimensional uniformity tests are either Chi-square tests or distance- 
based methods (Cressie, 1993). The first ones depend on the number and 
location of the quadrats (cells in which events are counted), whereas the 
last ones require numerous simulations. More recently, there has been some 
interest in e.d.f.-based methods and extensions of the Cramer-Von Mises 
test (Zimmerman, 1993) and the Kolmogorov-Smirnov test (Justel & al, 1997) 
to the [0, l] 2 case have been established. 

On the other hand, spacings theory, so useful for testing uniformity on R, 
remains almost unworked in higher dimensions even if one may think, as 
Zimmerman (1993) does, that distances from events to their nearest neigh- 
bours can be viewed as two-dimensional analogues of spacings. We shall just 
mention the results of Deheuvels (1983b) and Janson (1987) concerning the 
asymptotic distribution of the maximal multidimensional spacing, i.e. the 
volume of the largest square (or ball) contained in [0, l] 2 and avoiding every 
point of the pattern. 

A first application of spacings theory to CSR testing would be to test both 
x- and y-coordinates' uniformity using a spacings-based method. The rejec- 
tion of either leads to refuse the two-dimensional uniformity hypothesis. But 
we can never accept it as a bivariate distribution with uniform marginals 
need not be uniform. This makes necessary to take into account the joint 
distribution of the x- and y-coordinates. 

In this paper, following this idea, I introduce a new notion of two-dimensional 
spacings which is related to spacings based on x- and y-coordinates. This 
relationship then allows me to derive the limiting distribution, under the 
uniformity hypothesis, of a wide family of statistics based on these spacings. 
This is done by a direct decomposition method similar to the one Beirlant & 
al (1991) used for one-dimensional spacings. 

An application of this asymptotic result is developed by Cucala & Thomas- 
Agnan (2005). Two of these statistics, the variance and the absolute mean 
deviation of the two-dimensional spacings, are selected and used in practice 
to test for CSR. A multiple procedure is adopted to generalize the tests to 
point processes in any domain (not necessarily rectangular). Then the power 
of these spacings-based tests is compared to the power of existing tests us- 
ing real and simulated data sets : they appear to be inferior for detecting 
regularity or clustering but more powerful for detecting certain types of het- 
erogeneity. 
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Fig. 1 - Two-dimensional spacings 

2 Spacings in [0, l] 2 

2.1 Definition 

Let U = ((£/?, Uf), ■ ■ ■ , (E/£_ ls U y n _S) be a point pattern in [0, l] 2 . 
Let U x = (I/f, • • • , E£_ x ) and = (C/f , • • • , C^_ x ). 

UZs < • • • < f^(n-i) are the order statistics corresponding to U x . 
C/L < • • • < U? n ^u are the order statistics corresponding to U y . 

Set US = Ul = C/f 0) = C/fo = and E£ = E# = ^ = Ufa = 1. 
One may define the spacings related to the pattern 

the x-spacings Df = Uf^ — U%_ x y i — 1, • • • , n , 
the y-spacings D y = U y j} - U y {j _ iy j = 1, ■ ■ ■ , n . 

A way to take account of how x- and y-spacings vary jointly is then to define 
the two-dimensional spacings as the areas formed by the grid in Figure ^ 

Vi 6 {1, • • • , n},Wj G {1, • ■ ■ , n}, A tJ = D*D y . 

2.2 Uniformity hypothesis 

From now on, I will assume the uniformity hypothesis H : the point 
pattern U is uniformly distributed in [0, l] 2 . 

Let (Di, ■ ■ ■ ,D n ) be the spacings corresponding to a (n-1) uniform sample 
on [0, 1]. Then it is easy to see that 
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r C(Df, ••-,£>*) =C{D U ■■■,D n ), 
{ C(Df,... ,Dv) = C(D u --. ,D n ), 
( (Df,-.. ,D*)AL(Df,... ,Dy). 



3 Asymptotic normality of additive functions 
of spacings 

3.1 Main result 

Many statistics based on one- dimensional spacings are additive functions 

n 

V n = Y,9{nD i ) 

i=i 

for a measurable function g. 

For the two-dimensional case consider 

n n 
i=l j=l 

The asymptotic normality of V n was proved by Beirlant & al. (1991) using 
the following distributional equivalence (Moran, 1947) 

(n^,---,^) -(§,•••,§) (1) 

where (E 1 , • • • , E n ) are independent exponentially distributed r.v.'s with mean 1 

1 n 

and E= - V E t . 

n 

i=i 

(2) 

I will use the same technique to prove the asymptotic normality of Vn under 
H . 

Introduce (X±, • • • , X n ) and (Yi, • • • , Y n ) two samples of independent expo- 

1 n 

nentially distributed r.v.'s with mean 1. Denote X = — } X; and 

n — ' 



n 
1=1 



1 

Y = — > y,. Then define the statistic 

3=1 



i=i j=i 
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From (jTjl . 14 and G n have same distribution. 

From now on I will assume g satisfies the following, where ip and ip are 
measurable functions 

g continuous on R + , (2) 
Eg 2 {X 1 Y 1 ) < oo, (3) 
Vt G R + , Bip : R -> R, Vn G N, Etp(XY) < oo 

and Vi < t ,Vx G R + , |#(ta;)| < <p(x), (4) 
Vt G R + , 3^ : R -> R, Vn G M, E^(XF) < oo 

and V/ - /„.¥./• ■. ,/(/ '' ) 

Denote 



g(t) 



< ip(x). (5) 



/x = E[g(X 1 Y 1 )}, 

r, = Cov(g{X 1 Y 1 ),g(X 1 Y 2 ))=Cov(g(X 1 Y 1 ),g{X 2 Y 1 )), 
c = Cov(g(X 1 Y 1 ),X 1 ) =Cov(g(X 1 Y 1 ),Y 1 ). 

To justify my decomposition, I make the following argument. 

Using the same Taylor-expansion as Proschan & Pyke (1964), it appears that 



n 



-3/2 



(G n — r?\i) could be asymptotically equivalent in distribution to 



n n 



3/2 E E^Wi) ~ » ~ ~ ^s'Wj) WS) - (y - lypf^xx^)} 

i=l j=l 

n n _ n n 

E E {^(^ y i) -m-(^-i+^-i)(^EE wi)) }. 



n 

i=l j=l 
n n 

-3/2 



n 

i=l j=l fc=l !=1 

By partial integration it follows that 

-. n n 

Eg'(X l Y 1 )(X 1 Y 1 ) = c =s> — V V^^)^) c. 

7T, Z z — ' n— >oo 
fc=l £=1 

That is why it seems useful to decompose G n as follows 
■(G n -n u) = — - + 



3/2 VT" ^/ n 3/2 1 n 3/2 

where S n = ^ E^^i) " ^ " C ( X * " X ) " c (^' ~ ^ 



i=l j=l 

' Xi Ya 



and i?„ = EE KTF) ~ ^^)] + c ^ ~ ^ + cr ^ F " ^ 



i=l j=l 



s 

As —5-^ is a two-sample U-statistic with mean and limiting variance 
o 2 = 2(j] — c 2 ), one has from Van Der Vaart (1998) 



S n d 



n 3/2 

I will prove in the next section that 



AA(0, 



12° n^oo 



which will yield the following result. 

Theorem 1 Assume g satisfies J3J), (0) and J3J). Then 

1 



TV 



3/2 



{G n - n 2 /i} -=-> AA(0, a 2 



3.2 Behaviour of the remainder term 

Using the in< 
and Y, one gets 



X- Y- 
Using the independence of and X, as well as the independence of 



^C^»l) T-l I r) rji . rp 

5 — J l,n + z 1 2,n i 1 3,n 



rr 



where T\ n = — 
n 



E M§ f)} " 2E M# F> (XlFl) ) + E ^ X ^ 



x y/J r vx y 

// - I Lf /.V, V;\ f X x Y 2 



(6) 



7 - ^HHtfJHtfJ) - 2E 14tfJ^ x ^} 



n 



and T 3;n 



n 



X\Y\\ / X 2 Y 2 



-2c 2 . (8) 



+E^(X 1 y 1 )r 7 (X 1 y 2 ) ,(7) 



^y 



3.2.1 Preliminary results 

The marginal and bivariate densities of the spacings Di are given by Pyke 
(1965) and lead to 



9 fX 1 Y 1 \ (n-1) 2 f n f n „ . f i\n-2/ y\™- 2 , , 
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fX 1 Y 1 \ (MY2\\ 



X1Y2 

X Y 



(n- l) 2 (n-2) 



n rn rn—y 




g(xy)g(xv) 



>0 JO JO 

(* - lY X 1 - ^^y'dvdydx, (10) 



2 rn rn rn—x rn—y 



^mxfi)} - ( - i) :r 2)i / / / / 9 »» 




(l - ^) n "(l - ^y-'dvdudydx. (11) 



n-3 



x, t 



Using the independence of -=r and X, one finds 



E 



[n-lf rn rn 



rr 



Jo 



g(xy)E{g(xXyY)} 



Xi Yi Y^ll 
X' Y' Y )\ 



(n- l) 2 (n-2) 



rr 



n /•« rn—y 




g(xy)E{g(xXvY)} 



x 



n 



n-2 



JO JO 

y + -u^™- 3 
n 



dv dy dx, 



(12) 



(n- l) 2 (rr-2) 2 /•"-* /• n -«' 



X\ X 2 Y 1 Y 2 



rr 




Jo Jo 



' y' y /J 
y(x J/ )E{^(«x u y)} 







(l - ^) n " 3 (l - ^vdudydx. (13) 



n / v n 
The following lemma is also needed. 

Lemma 1 If g continuous on M + * and Eg 2 (X 1 Y 1 ) < 00, then \/t e [0, +oo[ 

lim E#(tX?) =g(t). 
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Proof : 

Denote by f n (u,t) = f^y-^ e ~ W * the 

common density of tX and tY. 

poo poo 

Eg(tXY) = Eg(\/ixVtY) = / g(uv)f n (u,\/i)f n (v,\/i)dudv 

Jo Jo 

g(uv)f n (u,y/i)f n (v,y/i)dudv + / / g(uv)f n (u,y/i)f n (v,y/i)dudv 

D J Jd 

where D = {{u, v) G R +2 ; y/i/2 <u< 3y/i/2; y/i/2 < v < 3y/i/2}. 

It is easy to see that : 

\/tX — ^ y/i and \ftY — ^ y/i. 

n^oo n— »oo 

Introduce the function y> : D —>■ R 

-> p(ar,y) = flf(:cy). 

From (0) the function <£> is continuous and bounded on D , so by the Helly- 
Bray theorem one concludes 




g{uv)f n {u, Vi)f n (v, Vi)dudv = Eip(VtX, VtY) ► <p(Vi, y/i) = g(t). 



D 



It remains to prove that f f n g(uv)f n (u, y/i)f n (v, \/i)dudv ► 0. 

K and K are two constants. From Beirlant & al (1991), one has Mn > 1 



n— 1 



f n (u,Vt) < Kt-Wn^^e 1 -** 

SO : (it, V) G D f n ( U , Vt)f n (v, y/i) > 0. 

n— >oo 

Moreover : n > n > 5 f n (u, \/i)f n (v, y/i) < Kf no (u, \/i)f m (v, y/i). 
And taking m > 16, we get, Wu G [0, y/i/ 2] U [3y/i/2, +oo] 

i-«M) m <e -«M 

=>• / no (w, < K^ 2 r x/2 ri /2 e~ u if we take n > 16^ + 1. 
As, from ©, («, v) — > g(uv)e~ u e~ v G £ 1 (M 2 ), one gets 

g(uv)f no (u,Vi)f no (v,Vi)dudv < oo. 

Lebesgue's dominated-convergence theorem leads to the conclusion. ■ 
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3.2.2 Behaviour of T l n 

V(x,y) G [0,n] 2 , 

|^ 2 (xy)(l-x/n) n - 2 (l-y/n) n - 2 | < ^ 2 (xy)e 4 e- a: e- 2/ G /^([O,™] 2 ) from ©. 
So applying Lebesgue's dominated-convergence theorem to (jHJ) leads to 

/ X Y \ Z" 00 Z" 00 
limEj 2 ^ = / ^(^)e- a, e-»dyd a : = E(/ 2 (X 1 y 1 ). (14) 

n^cxD v a y / y Jo 



By Cauchy-Schwarz inequality one gets 

= Eg 2 (X 1 Y 1 ) from (JUD 

< oo. (15) 

©, (HI) and (HSJ) lead to 

lim Ti n = 0. 



3.2.3 Behaviour of T 2 n 

V(x,7/) G [0,7i] 2 , We [0,ra- y], 

\g(xy)g(xv)(l—x/n) n ~ 2 (l — (y+v)/n) n ~ 3 \ < \g(xy)g(xv)e 5 e~ x e~ y e~ v \ G 
£*([(), n] 2 x [0,n-2/]) from JHJ). 

So applying Lebesgue's dominated-convergence theorem to (|TU|) leads to 



Introduce the function h~ 



t -> /i n (t) = E#(tXY~). 



Lemma 1 gives 

g(xy)h n (xv)(l - x/n) n_2 (l - (y + 7;)/7i) n ~ 3 > g(xy)g(xv)e~ x e~ y e~ v . 



Denote t\ G K + *. From (jSJ), one gets 
37/>:M^R,W>ti,Va;GM + 



0$ 



< ?/>(x) 



Kit) 




g(t) 






'0 

POO ^oo 



/ n (x, 1)/„(t/, l)dydx 



< 




i/>(xy)f n (x, l)f n {y, l)dydx = E^XY ) = *(n). 



o jo 



So, V(x,t>) G R +2 ,xv > ti =>- \g(xy)h n (xv)(l-x/n) n ~ 2 (l-(y+v)/n) n 3 | < 
\g(xy)g(xv)^(n)e 5 e~ x e~ y e~ v \ G /^([O,™] 2 x [0,n — y]). 

From (jlj), one gets 

3yj : R -> R,Vt < ti,Vx G R + , \g(tx)\ < <p(x) 
V* < *i, |M*)I < Ey?(X?) = $(n). 

So, V(x,v) G M +2 ,xf < ii \g(xy)h n (xv )(l-x/n) n ~ 2 (l—(y+v)/n) n ~ 3 \ < 
\g(xy)$(n)e 5 e- x e- y e- v \ G /^([O,™] 2 x [0,n-y]). 

So applying Lebesgue's dominated-convergence theorem to ()12jl leads to 

jim e{^(^^)^(x 1 f 2 )} = E 5 (x 1 y 1 ) 5 (x 1 y 2 ). (17) 

(JZD, (dnj) and (H7|) lead to 

lim To .„ = 0. 



-2,n 
ri —>oo 



3.2.4 Behaviour of T 3 , n 
From 1)11)1 ■ one gets 
(n-l) 2 m r /XxFa /X 2 F 2 



to 

( n _l)4( n _ 2 ) 2 p Z" 00 p p . . . ./ 3(a; + M )-(x + n) 2 /2 
(l + % + V) ~ n (y + v)2/2 )e-*e-ye-e-"dvdudydx + I n 



to j rr 
- 2 ^"-') 4 ("-^) 2 E[XiWl)1 _ 2 (^-l)V-2)^ 2[Xig(Xiyi)] + ^ 

where J„ = 

(n- l) 4 (n-2) 2 r P f n [ n ~ x [ n ~ y / x + «\™-3/ y + t;x«-3 



n 5 




'0 ^0 ^0 ^0 

~oo /.oo /.oo . 3(x + to) - (x + to) 2 /2 



dvdudydx — I \ \ \ g(xy)g(uv)(l + 
-'0 </o io 



n 



■ 1 + 3fa + ,)-fa + ,)V2x e _„ e _, 
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From (JT3J), one gets 

( n — 1 ^{n — 2) 2 / '™ / ' n /' n_:r , /" :x /" :x 

n 5 





g(ab)f n (a,u)f n (b,v)da db 



'0 JO JO JO X J0 Jo 



y(a;y) ^1 j ^1 — - — dvdudydx 



n 5 




o Jo Jo Jo 



U n y / y _|_ y \ n-3 
/ n (6,u)^l J dvjdbdadydx. 

Now the first integral in braces is equal to 



n — 1 V n) \ (n — x)(n — 2) 



(n - l) 2 (n - 2) 2 [ n f n f°° f°° , , , / s\«-3 



"" Jo Jo Jo Jo ra 

ra— 3 



e 



■»/<-*) ( 1 + ™ ) e -> (i _ £ Y 1 e -^/(-y) 

V (n — x)(n — 2) / V nJ 



(n — x)(n — 2) 

1 + r ) db da dy dx 

(n-y){n -2) J 

( n _l)2 (n _ 2 )2 fn fn foo foo 3 

' ' ' ' g(xy)g(ab)e 1 



n3 J"o J"o J"o J"o v /( 



n/ \ n — x (n — x)(n — 2) 

by Tib \ 

1 h r I db da dy dx + J n 

n-y {n-y)[n-2)J 

where J n = 

( n _l)2 (n _2)2 fn fn foo foo 3 yV -3 



n 3 




Jo Jo Jo 



g{xy)g{ab)e-\l-l) n ~ ^(l - £)' 



(; ax/(n-x) e -fey/(n-y) / j_ + nQ ^ A + ™ & 



(n — x)(n — 2) J \ (n — — 2) 



arc na \ / by nb , . 

1 h t r- r 1 h t -7 7 >dbdady dx. 

n — x (n — x) [n — 2) J \ n — y [n — y) [n — 2) / 
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n , . . ax ( x\ - 1 ax TT 
Substituting — II ) to , one gets that U n equals 



n 



n 



n — x 



("-l)>-^ + 6(n-l)'(n-2)' 



n 1 



2{n-l)\n-2f 




xg(xy)(l - 



+ »-"?»-V #a 



'0 ./0 

"n rn 



X 



n 



n— 4 



n 



n— 3 



dx 



rr' 




o 



g(xy)(l--) (l - -j dydx 



n 



(n- l) 2 (n-2)< 



■fiE[X^g(X 1 Y 1 )] +J n + K n + Oin- 1 ) 



(19) 



where K„ 



(n-l) 2 (n-2f 



fOO POO 



-ft 




JO 



g(xy)[l 



/• \ " 3 f y\ "-■'> 



n 



n 




g(xy)e- x e-y(l + 



3x - x 2 /2 



n 



3y-y 2 /2 



'o jo 

So, from (fSf. IjIBjl et (fT5j) we get 

T 3 ,„ = A n + 0(n- 1 ) + I n + J n + K n - 2c 2 



n 



dy dx 



dy dx 



(20) 



where A n = (n - 8)/i 2 + 12//E[Xi#(Xin)] - 2 / uE[X 1 2 r/(X 1 y 1 )] - 2E 2 [X 1 r/(X 1 y 1 )] 
- 2(n - Q)/j 2 - \2ijE\XxgiXxYt)} + 2 / uE[X 1 2 c/(X 1 y 1 )] + (n - 2)/j 2 



+ 4E[X l5 (X 1 Y 1 )] 
- 4E[X 1 < 7 (X 1 F 1 )] 




xg(xy)[l 



'0 JO 

en rn 



r \ " ; / ^ //\ 



dx 




o ./o 



x \ rt ~ 4 / . ?/ \ n ~ 3 



#(xy) 1-- 1-- dydx + 0{n~ l 



n 



A r 



2/i 2 - 4 M E[X 1 #(X 1 y 1 )] + 2E 2 [X 1 c/(X 1 r 1 )] 



2c 2 . 



(21) 



It now suffices to show that I n + J n + K n = o(l). 
A Taylor-expansion leads to 



Vx e [1, y/n\ : 1 - - 
* n 



n— 3 



1 + 3X x2/2 +0(n- 2 (x 2 + x 4 ))l. (22) 



n 



Denote b n = 41ogn. Choosing n large enough gives 



Vx E [b n ,n], (1 



x 



n 



n— 3 



< e" 



(23) 
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By Cauchy-Schwarz inequality one gets 

,3 



n 



[n - l) 2 (n -2) 2 /i 



\K n \ 




g(xy) 



Jo 



nJ \ n, 



- e~ x e- y (l + 



3x-i 2 /2\/ 3y-y 2 /2 



n 



n 



dy dx 



< \Eg 2 (X 1 Y 1 )] 1 ^ 




- e~ x e- y 1 



'o </o 

.,2 



nJ \ ri/ 



3s-x 2 /2w 3y-y 2 /2 



n 



n 



e x e y dy dx 



1/2 



where x + = x if x > 0, elsewhere. 
Denote _E n as the double integral 



oo poo 




JO 



1 + 



1 - 



X 



n 



n 



n—3 



e~ x e- y 



3x-x 2 /2\( 3y-y 2 /2\i 2 



n 



71 



e x e y dy dx. 



Using (|2*2*j) and (j2Sj) and following the technique of Does & Klaassen (1984), 
one can prove : E n = 0(r?T 4 ) 

K n = 0(n- 1 ). (24) 



The same arguments are used for 

I n = 0{n- 1 ) 



(25) 



Using the inequalities 1 — z > e z > 1 — z + z 2 /2, z > 0, the expression in 
braces in the definition of J n is bounded below by 



a 2 nx 



n 2 a 2 bx 



(n — x) 2 {n — 2) (n — 2) 2 {n — x) 2 {n — y) 
b 2 ny n 2 b 2 ay 



(n — y) 2 (n — 2) (n — 2) 2 (n — y) 2 (n — x) 
and bounded above by 
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a 2 x 2 na 2 bx 2 
+ 



+ 



2(n — x) 2 2(n — 2)(n — x) 2 (n — y) 
a 2 b 2 x 2 y 2 na 2 b 3 x 2 y 2 



4(n — x) 2 (n — y) 2 A{n — 2){n — x) 2 {n — y) 3 
na 3 x 2 n 2 a 3 bx 2 
+ 2(n-2)(n-xf + 2(n - 2) 2 (n - x) 3 (n - y) 
na 3 b 2 x 2 y 2 n 2 a 3 b 3 x 2 y 2 



+ 
+ 



4(n — 2)(n — x) 3 (n — y) 2 A{n — 2) 2 {n — x) 3 (n — y) 3 

b 2 y 2 nb 2 ay 2 
2(n — y) 2 2(n — 2)(n — y) 2 (n — x) 
b 2 a 2 y 2 x 2 nb 2 a 3 y 2 x 2 



A{n — y) 2 (n — x) 2 A{n — 2){n — y) 2 (n — x) 



3 



nb 3 y 2 n 2 b 3 ay 2 



2{n-2){n-yf 2(n - 2) 2 {n - yf(n - x) 
nb 3 a 2 y 2 x 2 n 2 b 3 a 3 v 2 x 2 

+ -T, =77 ^TT7 7T + 



4(n - 2){n - yf(n - x) 2 4(n - 2) 2 (n - y) 3 {n - xf 
Hence one gets 

3 pn pn poo poo 



J yl J dbdadydx 



+ 



n rn poo poo 




n(n 2) 2 J Jq J J 

J J dbdadydx 



\g{xy)g(ab)\ab 2 ye 



=> J n = Gin' 1 ) 
(PU|1. (EU), (El, (ESD and (j2SD lead to 

lim T 3i „ = 0. 

71— >QO 



14 



References 



Beirlant, J., Janssen, P. and Veraverbeke, N. (1991). On the asymp- 
totic normality of functions of uniform spacings. The Canadian Journal 
of Statistics, 19, 93-101. 

Cressie, N.A.C. (1993). Statistics for Spatial Data. Wiley, New York. 

Cucala, L. and Thomas- Agnan, C. (2005). Tests for spatial randomness 
based on spacings. Submitted to The Canadian Journal of Statistics. 

Deheuvels, P. (1983a). Spacings and applications. Probability and Sta- 
tistical Decision Theory. Volume A (F. Konecny, J. Mogyorodi and W. 
Wertz, editors). Reidel, Dordrecht, 1-30. 

Deheuvels, P. (1983b). Strong bounds for multidimensional spacings. 
Zeitschrift fur Wahrscheinlichkeitstheorie und Verwandte Gebiete, 64, 
411-424. 

Does, R.J.M.M. and Klaassen, C.A.J. (1984). The Berry-Esseen theorem 
for functions of uniform spacings. Zeitschrift fur Wahrscheinlichkeits- 
theorie und Verwandte Gebiete, 65, 461-472. 

J ANSON, S. (1987). Maximal spacings in several dimensions. Annals of 
Probability, 15, 274-280. 

Justel, A., Pena, D. and Zamar, R. (1997). A multivariate Kolmogorov- 
Smirnov test of goodness of fit. Statistics and Probability Letters, 35, 
251-259. 

Le Cam, L. (1958). Un theoreme sur la division d'un intervalle par des 
points pris au hasard. Publications de ITnstitut de Statistique de I'U- 
niversite de Paris, 7, 7-16. 

Moller, J. and Waagepetersen, R.P. (2004). Statistical Inference and 
Simulation for Spatial Point Processes. Chapman & Hall/CRC, Boca 
Raton, Florida. 

Moran, P.A.P. (1947). The random division of an interval. Journal of the 

Royal Statistical Society Series B, 9, 92-98. 
Proschan, F. and Pyke, R. (1964). Asymptotic normality of certain test 

statistics of exponentiality. Biometrika, 51, 253-256. 

Pyke, R. (1965). Spacings. Journal of the Royal Statistical Society Series 
B, 27, 395-449. 

Van Der Vaart, A.W. (1998). Asymptotic Statistics. Cambridge Univer- 
sity Press. 

Zimmerman, D.L. (1993). A Bivariate Cramer- Von Mises Type of Test for 
Spatial Randomness. Applied Statistics, 42, 43-54. 



15 



